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Abstract 

We present a numerical study of type IIB supergravity solutions with varying 
Ramond-Ramond flux. We construct solutions that have a regular horizon and 
contain nontrivial five- and three-form fluxes. These solutions are holographically 
dual to the deconfined phase of confining field theories at finite temperature. As 
a calibration of the numerical method we first numerically reproduce various an- 
alytically known solutions including singular and regular nonextremal D3 branes, 
the Klebanov-Tseytlin solution and its singular nonextremal generalization. The 
horizon of the solutions we construct is of the precise form of nonextremal D3 
branes. In the asymptotic region far away from the horizon we observe a loga- 
rithmic behavior similar to that of the Klebanov-Tseytlin solution. 



1 Introduction 



The AdS/CFT correspondence provides an alternative way to study the strongly 
coupled regime of gauge theories via gravity duals. The original statement of the 
AdS / CFT correspondence identifies M = 4 supersymmetric Yang-Mills with IIB strings 
on AdS$ x S 5 PUI2|- At finite temperature the gravity side is described by nonextremal 
D3 branes and the qualitative matching of the properties was one of the key observa- 
tions in the understanding and eventual formulation of the AdS/CFT correspondence 
jSj- Other interesting aspects of finite temperature theories, as seen by the AdS/CFT 
correspondence were discussed by Witten In particular, the Hawking-Page phase 
transition in the gravity side was related to the confinement /deconfinement transition 
on the field theory side. 

In a series of papers, Klebanov and collaborators [21101111013110] developed a gravity 
theory that encodes interesting field theory phenomena like chiral symmetry breaking 
and confinement. The class of Klebanov-Tseytlin (KT) solutions are solutions of IIB 
supergravity with nontrivial metric and F 3 ,H 3 and F 5 forms. A phenomeno logically 
very attractive property of the supergravity solution is that it encodes the logarithmic 
running of a gauge coupling in field theory. It does so via a varying E>2 field which is 
compensated by a constant F 3 flux through a 3-cycle. The five-form, which is constant 
in most solutions, varies according to the Bianchi identity dF 5 = H 3 AF 3 and generates 
a varying flux as F 5 depends on the radial coordinate. The supergravity solution 
therefore has varying flux. 

The finite temperature phase of such theory is certainly very interesting and has 
been tackled in various papers including fTJ El Ej In particular, reference (El 
constructed a perturbative solution whose regime of validity is restricted to high tem- 
perature and small value of the F 3 flux: F 3 = P. Knowing the solution only 
asymptotically in the radial coordinate and for a specific regime of parameters pre- 
vents us from extracting the full thermodynamics and from being able to understand 
possible phase transitions. In this paper we construct numerical solutions with regular 
finite- area horizons and nonvanishing values of the five and three-form fluxes. In gen- 
eral we find that the solutions have horizons of the form of the nonextremal D3 branes 
and in the other asymptotic region are particular generalizations of the KT class of 
solutions with the expected logarithmic behavior of the fields. 

Although we have motivated the study of these backgrounds from the AdS/CFT 
point of view it is worth mentioning that from the purely gravitational point of view 
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these are novel backgrounds. Black holes with varying flux have not been studied 
outside the context of AdS/CFT. These black holes pose interesting problems in the 
gravitational sense. For example, the definition of conserved charges is ambiguous [T3] . 
Recently, a better understanding of their thermodynamics and other properties was 
provided in f5J E>] . Our study might also be of interest to numerical relativists given 
that the qualitative properties of the solution strongly depend on the parameters and 
initial conditions. For instance, we intuitively expect the existence of the horizon to 
depend on various parameters including the temperature itself. For very low tempera- 
ture we expect that the solution will either develop singularities or be suppressed with 
respect to the Klebanov-Strassler background. Numerical methods are particularly 
suited to tackled these kind of problems. 

Strictly speaking, our result is relevant for a finite temperature field theory that 
has jV = 1 supersymmetry at zero temperature. Interestingly, this theory is confining 
and we expect that this line of research might eventually provide a framework for 
understanding RHIC physics. 

The organization of the paper is as follows. In section El we present the Ansatz for 
the metric and form fields involved in the solution, we also reduce the ten-dimensional 
problem to a system of equations depending only on the radial coordinate. Section 
El also contains the analytic form of various known solutions. In section El we present 
a description of the numerical method and study numerically the known solutions 
presented in section |21 The error of the numerical approach is estimated as the square 
root of the square difference between the known analytical solutions and our numerical 
solutions. We establish that typically this error is \ < 10~ 13 . Section 0] contains 
our main result which is a numerical solution with regular horizon of finite area and 
three-form flux turned on top of the standard five-form. We conclude in section In 
appendix El we discussed the geodesic behavior of various typical solutions to show 
how these features can be used to identify properties of the new solution. 

2 Review of gravity backgrounds with varying flux 

In this section we describe the setup from the ten-dimensional point of view and reduce 
it to a system of ordinary differential equations. The work presented in this section 
is largely a review of [T2*l IT3j and we refer the reader to those papers for the original 
presentation and details. Our intention is to be self-contained since we are going to 
use these analytical solutions to calibrate the numerical method. 



2 



The main idea in constructing the solution is to replace S 5 by a five dimensional 
manifold known as T 1 ' 1 which is topologically a product of a 2- and a 3-sphere. This 
manifold will be parametrized by coordinates (■0, 9i, 0i, 62, $2)- The Ansatz for a 
nonzero temperature generalization of the KT solution contains various fields. The 
construction of the Ansatz follows directly the one presented originally in ^2] • For the 
metric we consider a generalization consistent with the U(l) symmetry generated by 
-0-rotations. 1 The Ansatz in question depends on four functions x, y, z, w of the radial 
coordinate denoted by u: 

ds 2 = e 2z {-e~ 6x dX 2 + e 2x dX,dX l ) + e~ 2z dsl , (2.1) 

where 

dsl = e 10v du 2 + e 2y (dM 5 ) 2 , 
(dM 5 ) 2 = e-^el + e^iel+el+el+eD^e^dsl (2.2) 

The Funfbein is: 

1 11 

e i> — 9(^0 + cos9\d(f)i + cos02d(j)2) , te i = —j=d9i , = —=sm9id(j)i . (2.3) 
3 V6 v6 

The coordinate X° represents time and X 1 are the 3 longitudinal directions of the 
3-brane. The qualitative meaning of the metric functions x, y, w and z can be clarified 
from the Ansatz. The function x breaks the Poincare invariance in the four plane 
defined by (A , A*) and therefore describes the nonextremality. The function z multi- 
plies the four-dimensional and the six- dimensional subspaces of the metric by different 
factors and can be clearly interpreted as the warp factor. The function y basically 
amounts to a choice of the parametrization of the radial coordinate, it introduces a 
natural ambiguity that we discuss later on. Finally, the function w describes how the 
U(l) fiber is squashed with respect to the four- dimensional base. 

To reinforce the meaning of the metric functions x,y,z and w, it is convenient to 
use that the intuition for these gravity solutions has been developed in the context of 
D3 branes. We recall that the relation with D3 branes can be established as follows: 

ds 2 = h-^^i-g^dX 2 + dX t dX t ] + h l l 2 (p)\g 2 \p)dp 2 + p 2 dsl] , (2.4) 

the gauge theory [H] this symmetry is identified with the U(1)r. Restoring this symmetry at 
high temperature is understood as chiral symmetry restoration |13| . 
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with the identifications 

h= e - 4 *- 4x , p = e y +x+w ^ gi = e -8x^ e Wy+2x du 2 = g^{p)dp 2 . (2.5) 

Note that in the absence of nonextremality (x = 0), h = e Az which shows that z is 
truly the warp factor. In the absence of U(l) fiber squashing, that is w = 0, one has 
e 4y = p 4 = 4-, the transverse 6-d space is the standard conifold with M 5 = T 1,1 ; this 
shows that y amounts to a choice of the radial coordinate. Small values of u correspond 
to large distances in p and vice versa. 

The Ansatz for the p-form fields is dictated by symmetries and will be taken to be 
as in the original KT solution [Hj: 

F 3 = Pe^p A (ee 1 A — e# 2 Ae^) , 
B 2 = f(u)(e 6l A - ee 2 Ae^) , 

F 5 = T + *F , T = K(u)e^ A e 6l A e 01 A e e , 2 A e 02 . (2.6) 

Note that the form of F 3 is such that it describes a constant flux along a 3-cycle, that 
is, J S£j F 3 = P. In some other notation this is called the number of fractional D3 
branes. As briefly alluded to in the introduction, the Bianchi identity for the 5-form, 
d* F 5 = dF 5 = H 3 A F 3 , implies 

K{u) = Q + 2Pf(u) . (2.7) 

That is, in the presence of 3- form flux (P ^ 0), the flux of F5 varies with the radius. 
The fact that K(u) depends on the coordinate u is very novel and has interesting 
physical implications. 

One can summarize the presence of flux in the background by noticing that all 
p-form fields are completely determined by a constant P and a function of the radial 
coordinate f(u). 

2.1 The system 

The system of equations was obtained in ^2] by reducing the problem to a one dimen- 
sional effective action for the radial evolution. We refer the reader to ^2] for details of 
the derivation. The idea is to plug the metric and all the forms in an effective action 
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for IIB supergravity and integrate with respect to all variables except the radial coor- 
dinate u. One is thus lead to an effective classical mechanical system. The simplest 
equation is for the nonextremality function x: 



x' — , x = au , a = const. (2.8) 

The reason for such a simple equation is that, as explained in jT2], it does not appear 
in the effective one dimensional Lagrangian, except for its kinetic term. 

The other functions y, w, z, f and $ are to be determined from a coupled system of 
equations: 

lOy" - 8e Sy (6e- 2w - e - 12w ) + $" = 0, 
lOw" - 12e 8y (e~ 2w - e~ 12w ) - $" = 0, 

$" _|_ e -*+4z-42/-4jH/j/2 _ e 2$+8y+8«) p2^ _ g 
4 Z " - (Q + 2Pf) 2 e 8z - e -*+4«-4»-4««(/'2 + e 2<f+8 y +8« ) p2) = Q> 

The integration constants are subject to the zero-energy constrain T + 1/ = 0, i.e. 

5y '2 _ 2Z ' 2 - hw' 2 - i$' 2 - _ e 8y (6e -2» _ e -12 M) 

+ l e $ +4z+ 4 J/+ 4 U ,p2 + l e 8^g + 2 pjY _ 3a 2 = Q (2 _ 1Q) 

As mentioned while introducing the Ansatz, the function y amounts to a choice of the 
radial coordinate. This is relevant for understanding the dimensions of all quantities. 
Note that in particular, the system (J2.9|) and the constrain (|2.1U|) are invariant under 
e y — > L e y and u — > L^ A u if we assume that Q — > L$Q, P — > LqP, a — > L^a and 
/ — > LqJ. We therefore express all dimensionfull quantities in units where Lq — 1. 



2.2 Some analytic solutions 

In this section we review some known analytic solutions to the system. Our goal is to 
develop the necessary intuition into the properties of the solution we are seeking. The 
existence of exact analytical solutions provides us with the unique opportunity to test 
the numerical method. The numerical treatment of these solutions will be presented 
in section El and this section can be considered as preparatory. 
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2.2.1 The blown up conifold 

One of the simplest class of solutions to the above system has the form of 1R 3,1 x CY, 
where CY stands for a Calabi-Yau space and more concretely a Ricci-flat space. The 
simplest solution corresponds to the conifold and is defined by 

e * = i-, (2.11) 

and all other functions are zero or constant in the case of the dilaton. In this class of 
solutions we also find the blown up conifold 2 JTHj. This solution is usually written as 

dsl = K^dr 2 + r 2 [Ke\ + e\ + e% + e 2 2 + eJJ , (2.12) 

where 

K = 1 - ^ = e- 10w , r = e y+w , b < r < oo. (2.13) 

This space has two nontrivial functions: y and w. The direct relation between r and u 
is given by: 



r 5 (l - -b) 



In general the relationship between the coordinates u and r is nontrivial but for large 
r we recover small u ~ 1/r 4 . The meaning of w is clearly to squash the U(l) fiber 
with respect to the four- dimensional base. Note that as r — > oo the squashing vanishes 
w -> 0. 



2.2.2 Non-Extremal D3-brane solution: Singular and Regular 

The general form of the nonextremal solutions contains a singular horizon. We should 
consider it here to show how the singularity looks in the numerical analysis. Imposing 
regularity of the horizon of the solution leads us to the standard nonextremal D3 brane 
(henceforth we will call it just the 'D3 solution'). 

Let us consider the general system of equations (|2.9J) with P = and / = 0, that 
is, with no varying flux and also with $ = and w — 0. Then, we are left with ()2.8|) 
and the following system: 

y " - Ae 8y = , z" - ^Q 2 e Sz = , (2.15) 
2 An adjustment in the periodicity of tp is required. 
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i.e. 

x ' = a, y' 2 = b 2 + e 8y , z' 2 = c 2 + q 2 e 8z , q = , (2.16) 
with the integration constants a, b, c related by the zero-energy constraint 

5b 2 - 3a 2 - 2c 2 = . (2.17) 

Assuming that a, b, c > to preserve asymptotic conditions, the general solution is: 

e*y = — b —- , e 42 = . UA \ , e ix = e Aau , (2.18) 

smn Abu q smn Ac[u + k) 

where k is defined by e 4cfc = o _1 (\/g 2 + c 2 + c) = 7 . Note that the integration constant 
k implies just a shift of variable u. To meet AdS asymptotic conditions we must choose 
k = 0. 

In terms of the familiar D3 brane Ansatz ()2.4j) : 

P 4 = e 4 ^ +4 - = \ _ e _ 8b , , £1 = <?2 = 2 = e" 8 -, (2.19) 

h = e _ 4 ,-4x = ± e A(c-a)u^ _ g-Scu) _ 

2c 

At small m (large p) we have 

<? = !-— + ..., p 4 = + ... . (2.21) 



The general solution with arbitrary b and c reduces to the standard extremal D3- 
brane background only if we set b and c proportional to a, satisfying the constrain 
()2.10|) . As we have emphasized several times, for arbitrary values of b and c the solution 
is singular. To see why that is the case, let us consider a singular case, discussed at 
length in [12], to better understand the role of parameters in the singularity property 
of the solution. The simplest special case is c = where z satisfies the 1-st order 
equation z' = —qe 4z . One can solve the system (j2.9|) to get 

e "-^L- u - e " 42 -^- »=^- < 2 ' 22 > 

and thus 

p 4 = e 4y+4x = °_ e _ 8bu , h = e~ 4z ~ 4x = Aque- 4au , g = e- 8au . (2.23) 
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Note that p is not well-defined as a radial coordinate for b ^ a. This solution has a 
singular horizon at u = oo as described in [T2|. Using the expression for R mn kiR mnkl one 
concludes that the metric is singular. Most importantly, the area of the horizon, which 
is proportional to exp(— 2z + 3x + by) vanishes as u — > oo. Namely, it is proportional 
to w 1 / 2 exp(— (y/TE — S)u) — ► as u — > oo. One 0/ t/ie main criterium for the solution 
we are seeking is thus, to have a finite nonzero area of the horizon. We have seen that 
the choice of the constants a, b and c plays an important role in fulfilling this condition. 

For small u (large distances) and in the limit a — > 0, we still get the standard 
asymptotic extremal D3-brane behavior 



P 4 = -L + ..., g = l-?$ + ..., h = \ + .... (2.24) 
4u p 4 p 4 



Note that in the large— p asymptotic the solution is like the standard D3. Most of the 
problems, as was pointed out, are localized at the horizon. 

The standard non-extremal D3-brane solution corresponds to the case when 

b = c = a, (2.25) 

This condition, together with the constrain (|2.17j) means that regularity picks a line in 
the two dimensional surface which is the space of solutions. We will see that something 
similar happens in the solution we construct in section 0] The solution of the system 
()2.9|) takes the form: 

p 4j/ _ a dz _ a Ax _ Aau Icy 

smh Aau q smh Aau 

Note that near the horizon (u — > 00) 

y = y,-au + ^e~ 8au + 0(e~ Wau ) , z = z* - au + ^e~ Sau + 0(e~ l6au ) , (2.27) 

1, 1 2a . 

y*= 7 ln2a, z, = - In — , 2.28 
4 4 q 

We can verified that this solution has finite area of the horizon. Indeed exp(— 2z + 

3a; + by) — > exp(— 2z* + 5y*) which is finite, g 1//2 (2a) 3//4 , as we approach the horizon 

(u — > 00). We see that the requirements on the asymptotics of the functions y and z 

in order to have a horizon with finite area are very stringent. 

A more recognizable form of this solution is given as 

ds 2 = h- 1,2 (gdXl + dXidXi) + h ll2 \g- l dp 2 + p 2 {dM 5 ) 2 } , (2.29) 
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p 4 ' ' 1 - e ~ 8au ' " p 4 

Let us clarify the relationship between the radial coordinate u and the more stan- 
dard coordinate p. Using 2a = p 4 , we have 

(2.31) 

This can be integrated to 

Note that in the domain of p < p < oo we have that u ranges in < u < oo. The 
position of the horizon which is finite in the p coordinate becomes infinite in the u 
coordinates. 

Some thermodynamics and universality of regular nonextremal D3 brane horizons 

For the standard nonextremal D3 brane the horizon is located at p = p , and the area 
is given by: 

A = Vu 5 R 2 pl, (2.33) 

where V is the volume due to the coordinates X 1 , R is the radius of AdS and appears 
in the harmonic function as h = R A / p 4 , u^, is the volume of T 1,1 . 

The natural temperature associated with the nonextremal D3 is obtained from the 
regularity of the Euclidean section: 

T=-%- (2-34) 



Of course, the local temperature is TJ oca / = TRpj a/ p 4 — p 4 and decreases as p — > oo. 
More explicitly, the nonextremal D3 brane metric is of the form 

ds 2 = ^ ((1 - f 4 )dr 2 + dXifUP) + y(l - ^ 1 dp 2 + R 2 du 2 5 . (2.35) 
Near the horizon one can introduce the following coordinate 

2 

P = Po(l + ^). (2.36) 

In terms of this radial coordinate in the limit of r — > the relevant part of the metric 
takes the form 

ds 2 = dr 2 + r 2 (j^drj . (2.37) 
9 



For the angular part to have period 2tt we obtain the above temperature as the inverse 
of the period: T = 1/(3. 

In the w-coordinates the area of surface defined by a horizon at u — constant is 

A = Vlu 5 exp (-2z + 3x + by) . (2.38) 

Given that the equation of motion for x has the general solution x = au we are forced 
into the following situation. If the horizon is at u — > oo, then in order for the area A 
to be finite we need the following asymptotics for z and y: 

z^aau + z*, y — > (3 au + y*, (2.39) 

with the condition that 

-2a + 5(3 = -3 (2.40) 

Note that the regular nonextremal D3 brane corresponds to a — (3 — — 1. The main 
claim is that: The existence of a regular horizon fixes the asymptotic behavior of the 
metric coordinates x,y and z. 

Similarly, one can obtain an expression for the temperature. Namely, the relevant 
part of the metric is 

ds 2 = e 2z - 6x dr 2 + e~ 2z+10y du 2 . (2.41) 
We introduce a new radial coordinate as: 

p = e z ~ 3x . (2.42) 



We can now rewrite the metric as: 



ds 2 = 



e 



-4z+l0y+6x 



\z' — 3x 



A2 



dp 2 + p 2 (e 2z - 5y - 3x (z' -3a) drf 



(2.43) 



Note that, again, the requirement of finite temperature fixes the large-u asymptotic of 
various metric functions to be 2z — 5y — 3x — * constant. Imposing absence of conical 
singularity we find that the temperature defined as the inverse of the period is 

rp = \ a ~ 3oj c 2z«-5^ (2.44) 

where we used that near the horizon the asymptotic form of z ~ aau + z* and similarly 
for y. 
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2.2.3 The Klebanov-Tseytlin background 

An important property of this solution is that it describes the asymptotic of most known 
solutions with P^O and for small values of u, see for example |H1 El H3 UH1 EH 120] • 
The KT solution can be represented as: 

x = , w = , $ = 0, 

e~ 4y = 4u, f = f -^\nu, (2.45) 

e~ 4z = K u - ^-u(\nu - 1) , K = Q + 2Pf , 

i.e. 

p2 p2 

e -^ = h = (Q + 2Pf + —)u-—u\nu. (2.46) 
For convenience, in the case of P ^ we can parameterize the function K as 

p2 

K(u) = -— ln(uL P 4 ). (2.47) 

As mentioned before the KT solution is an attractor of sorts. Let us clarify to what 
extend it is generic. The main claim is the following: Any solution of the class we 
considered with P^O will asymptote to the KT solution if we impose a constant value 
of the dilaton in the asymptotic region. We will not proceed to give a formal proof 3 
of the above statement, instead we will show how it comes about. If we impose the 
conditions of $ = and w = in the system of equations we find that the equation 
for y yields: 

y" = 4e Sy , (2.48) 
while the equation for the dilaton in ()2.9|) implies that 

f' 2 = P 2 e 8y . (2.49) 

One can explicitly solve this system and obtain that as u —>■ one has / ~ Inu. 

Let us present the ten-dimensional analysis which provides an interesting perspec- 
tive. For very large values of p we have that the metric becomes essentially a cone 
over T 1 ' 1 , in the effective one dimensional system language this is equivalent to having 
w = 0. The metric is thus 

dp 2 + p 2 ds 2 (M 5 ). (2.50) 
3 A rigorous analysis of fixed points of the system l|2.9|l will be presented elsewhere. 
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Imposing that we have P^O means that we are introducing a flux described by 

F 3 = Pn 3 . (2.51) 

If we also assume that the dilaton is constant we have and equation between the forms 
F 3 and H 3 which can schematically be written as I-F3I 2 = l-f^l 2 - This equation can 
also be understood as a consequence of having an imaginary self-dual G 3 = H 3 + iF 3 . 
Basically, it can be understood as a consequence of *6 F 3 — H 3 , where subscript means 
that we consider the Hodge dual in the six dimensional space. This is a very simple 
equation to analyze. In fact, one has 

* 6 F 3 = * 6 pn 3 = P^n 2 . (2.52) 

P 

Thus, we determine that B 2 = PQ 2 lnp, where Q 2 is some 2-form Hodge dual to fl 3 . 
This logarithm in B 2 propagates to a logarithm in the metric via the equation for F 5 
which is solved by 

K{p) = Q + 2P 2 In p/ Po . (2.53) 

Clearly, this solution has a repulsive singularity around p S i ng = po e~ Q ^ 2p2 which is 
locate at small values of p meaning large values of u (see appendix El for more details 
on the causal structure). 



2.2.4 A singular generalization of the KT Solution 

This solution is obtained by assuming that the dilaton is constant and w = in ()2.9|) 
and allowing a 7^ 0. The solution for y takes the form of 

sinh46w ' ^ 
with b = a a/3/5. Knowing y allows us to determine the other variables and one finds: 

P 

f = /* — — In tanh 2bu , 
e" 4z = ^ + ^(Li 2 (-e- 4b ")-Li 2 (e- 4b ")^) , (2.55) 

and we recall that K{u) = Q + 2Pf(u). This solution was originally obtained in jTTj 
as a nonextremal generalization of the KT solution. It was subsequently studied in 
[T2*] . where it was established that it has a horizon which coincides with a singularity 
at it = 00 and that it reduces to a singular nonextremal D3 brane for P — > 0. 
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2.2.5 Perturbative solution with varying flux around nonextremal D3 

In this section we review a solution presented in |1S 3 . This solution is obtained as a 
first order correction in P to the nonextremal D3 brane background. The main idea is 
to perform a perturbation around the regular nonextremal D3 brane (P = 0), in the 
regime where P 2 /K* <C 1 and it serves as a small parameter. 

It is therefore convenient to rescale the fields by appropriate powers of P 2 , setting 

K{u) = K* + 2P 2 F(u) , = P 2 (f){u) , w(u) = P 2 uj{u) , (2.56) 

and 

y^y + P 2 i, e~ iz -> e' 4z + P 2 ( , i.e. z -> z + P 2 r] , ( = -Ae~ iz r] + 0(P 2 ) , 

(2.57) 

where y,z represent the pure D3-brane solution ()2.2.2|) e~ Ay = a -1 sinh 4cm, e~ 4z = 
x 31 sinh4aw, and £ and ( or r] are corrections to it. The system ()2.9|) takes the following 
explicit form: 

lO£"-32Oe 8 ^ + 0" + O(P 2 ) = 0, 
lOu" - 120e 8y u - <p" + 0{P 2 ) = , 

(e 4 ^ 4y F')' - K,e 8z + 0(P 2 ) = 
Arj" — 8K 2 e Sz r] — AK*Fe 8z — e iz ~ iy (F' 2 + e 8y ) + 0(P 2 ) = (2.58) 

The constraint 12.101 becomes 

lOy'g -4z'r]' -~e iz - iy F ,2 -40e 8y Z+~e' iz+i y + K 2 e 8z r] + ^K,e 8z F + 0(P 2 ) = . (2.59) 

The solution takes the form 

K(u) = K*- ^-\n(l-e- 8au ) , (2.60) 

= 0* + -^Li 2 (e- 8a "), 



' v y ' r * 24^* 
The rest of the fields can be expressed in terms of a new radial coordinate of the form 

v = l- e' 8au . (2.61) 



13 



In this case we obtain: 

« = ^ + i 0^((- 2 « + ("- 2 )Ml-«))lB, + (,-2)Li 2W ), 

1 = ^^-(''"'Ml-fJ+Lijfi'))- (2-62) 

Similarly, one obtains a simple equation for u (see ^3] for details). 

Regime of validity, extrapolation and thermodynamics 

The above solution was constructed with the assumption that P 2 /K* <C 1. Since we 
considered only the linearized approximation there is an intrinsic limit on the values 
of u. Recall that for nonextremal D3 branes K{u) = Constant ~ which counts the 
number of D3 branes. For large values of u, the solution (|2.6U|) becomes 

K{u) ~ K* + ^-e~ 8au . (2.63) 

The second term is subleading for large u and thus in this region the solution remains 
valid. As we decrease the value of u, we reach a point where the first and second term 
in K{u) given by ()2.60|) are of the same order. This happens for 

u c = -— ln(l - e~ 2K *P 2 ) « L e -™*/P\ (2.64) 
8a 8a 

The value u c is small but nonzero. Altogether, the regime of validity is: 

P 2 /K, < 1, and u > u c = — e - 2K * /p2 . (2.65) 

8a 

Nevertheless, the authors of [TH] decided to explore the small u regime and found 
that 

p2 

K(u) ~K*-— ln(8aw). (2.66) 

Amazingly, this form of the solution resembles the KT solution. The main idea of (THj 
is to match the u — > asymptotics of the solution (j2.66J) with the KT solution (see 
equation 3.8 of [13J) 

K KT = ~H^%) (2.67) 
Matching these two solutions, that is ()2.66|) and ()2.67|) we find that (5.15) of |T3]) 

8aL" 4 = e 2K 'l p \ (2.68) 
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This is the relationship between and a. Finally one has the relation between a - the 
nonextremality parameter and the temperature 4 . Formula (5.16) in ^3] is nothing but 
the temperature of nonextremal D3 branes as a function of the radius of the horizon 
a 1//4 and the number of D3 branes K*. Thus, we quote 

T ~ a^K; 1 / 2 . (2.69) 

Plugging this relation into ()2.68j) we find that to leading approximation 

p2 rp 

K *~T ln A' (2 ' 70) 

which is quoted in (5.44) of This expression is crucial in understanding the 

thermodynamics of this class of solutions. One basically conjectures that the entropy 
per unit volume satisfies 

S P 4 / T x 2 



In- . (2.71) 



VT 3 L% V A 

A similar formula was recently discussed in the context of holographic renormalization 
in in 



Not all is well with extrapolating the perturbed nonextremal D3 brane solution 

This perturbative analysis is a very valuable tool to understanding the thermodynam- 
ics. However, it has several shortcomings. First, the thermodynamics requires the 
understanding of the region where the perturbation breaks down. Second, one can 
confirm that there are obstructions to extending the solutions analytically past its 
regime of validity. One can check that the constrain is not satisfied. In fact, expanding 
the constraint (J2.59j) for small u gives a divergent term already at first order in P 2 /K Jf : 

— — - -a 2 + -a 3 u + Oiu 2 ) 
2u 3 3 

3 Numerical analysis of known solutions 

To search for numerical solutions we rewrite ()2.9|) as a system of first order differential 
equations. To ensure that the numerical solutions automatically satisfy the Hamil- 
tonian constrain, the corresponding equations were modified to include explicitly the 
information given by expression (|2.10j) . The resulting system contains ten coupled fields 
4 This is an approximate relation that ignores some issues of asymptotics. 
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(2.72) 



{y(u), z(u),w(u), f(u),y(u), z(u),w(u), &(u), f(u)} related by the corresponding 
non-linear first order differential equations: 

dy _ dz dw d$ ~ df ~ 

t = y, -r = z ' ~i~ = w ^ ~r~ = $ ' ~r = Ji t 3 - 1 ) 

du du du du du 

dy ~2 ^ ~2 o ~2 ^ ^2 

— = 2y z — 2 w <P 

du y 5 20 

+ ^ e 8 ! /( 6e -2-_ e -12^ + JL ( g + 2 p /) 2 e 8,_6 a2) (32) 

^ = 5y 2 - 2z 2 - 5w 2 - i$ 2 
du 8 

+ | (Q + 2 Pf) 2 e Sz - e 8y {6e~ 2w - e~ 12w ) + L*+^+VM«>p2 _ 3^ (3 3) 
8 2 
9 4 9 9 1 -9 
~du = -2f + -i 2 + 2fi. 2 + -4 2 

+ y/" (e- 2 " - e- 121 ") - ^ « + 2i7) 2 e & + ^a 2 , (3.4) 

rf$ 9 9 9 1 ~ 9 

— = -20y 2 + 8z 2 + 20w 2 + -$ 2 
du 2 

+ 4e 8y (6e- 2w - e - 12w ) - - (Q + 2P/) 2 e 82 + 12a 2 , (3.5) 



dj_ 

du 



($ - 45 + Ay + Awj f + P (Q + 2Pf) e *+4*+%+4™ (3^ 



where P, Q and a are the three parameters that determine the behavior of the solutions. 
To solve it we use a combination of the methods implemented in the Maple software to 
find numerical solutions of ordinary differential equations 5 . For most calculations we 
have used the seventh-eight order continuous Runge-Kutta method which, thanks to 
its adaptive scheme, provides a great control upon the output accuracy. In those cases 
where the stiffness typical of singularities was present, we switched to the Livermore 
Stiff Ode solver. This also allowed us to get some information about the stability of 
the obtained solutions with respect to the truncation error. 

The option of yielding the output of the numerical computation as a list of proce- 
dures is, in general, not very economical for long and cyclical calculations. However, 
it proved to be convenient for the study of those quantities depending on the metric 
functions and on the dilaton such as those describing the thermodynamics. 

To use a numerical solver, one first needs to set up its arguments with values pro- 
viding outputs as accurate as required. We recall that in an adaptive scheme the 
discretization of the independent variable domain is automatically refined until a the- 



3 Our Maple sheets are available upon request. 
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Figure 1: Plot of the difference between analytical solutions and the corresponding 
numerical outputs: M 3,1 x CY (red), Klebanov-Tseytlin (black), standard nonextremal 
D3 (green) and singular nonextremal generalization of Klebanov-Tseytlin (blue). 

oretical measure of the error goes below some tolerance previously fixed by the user. 
So, to set the optimal tolerance we compared the known analytical solutions presented 
in the previous section with numerical outputs, using the measure given by, 



Here vf n (u) and v™ um {u) stand for the value at u of each of our ten variables as given, 
respectively, by the analytical and numerical solutions. Next, we must determine when 
this measure can be taken as negligible. We say that \ is negligible if adding its 
value to one of the parameters entering our system does not lead to a qualitative 
departure from the analytical solution corrresponding to the unchanged value of the 
given parameter. For instance, starting with a very low tolerance and a = we 
gradually incremented a until we found a qualitative departure from the KT solution. 
Then, we increased the tolerance looking for the maximal value that keeps the just 
described situation essentially unchanged. The same procedure was carried out with 
P and the standard non-extremal D3 solution. This way we found that for a tolerance 
of 10~ 14 , Xmax = 1CT 10 can be safely regarded to be negligible. 

In figure ^ typical errors are presented for the following solutions: KT (black), 
singular nonextremal generalization of KT (blue), D3 (green) and M 3 ' 1 x CY (red). 
The error measure is given by Eq. (|3.7|) and the plots show that, after tuning up the 



X = 



y/Y, \vr{u)-vr m {u)}\ i 



1- ■ -10. 



(3.7) 
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Figure 2: Plot of the numerical output of the constrain for four study cases: M 3 ' 1 x 
CY (red), Klebanov-Tseytlin (black), standard nonextremal D3 (green) and singular 
nonextremal generalization of Klebanov-Tseytlin (blue). 



numerical method, the calculations typically yield reliable results for a wide range of 
values of the independent variable u in the interval (10~ 3 ,20), and of the involved 
parameters. 

In the absence of a pattern solution to check for the accuracy of the computations, 
we fortunately have the Hamiltonian constrain. Since our system was explicitly con- 
structed to yield only solutions satisfying equation (|2.1(Jj) . the necessary fulfillment of 
this condition could be readily tested for every case we solved. As an example, we 
present in figure 121 the numerical results for the constrain corresponding to the four 
study cases mentioned above. 

3.1 Understanding the numerical solution 

Next, we move into understanding the output of the numerical computation. Our goal 
is to establish the existence of solutions with a nonsingular horizon and therefore we 
need to develop the appropriate criteria at u — > oo. At the same time we aim to make 
a definite statement about the behavior of the solutions for small values of u. 

It is important to note that since this is a nonlinear system of differential equations, 
the asymptotic form of the solutions is typically very sensitive to the values of the ten 
'initial' conditions. So, to search for a solution with given behavior at small and at 
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large values of u by just guessing the conditions at any of these boundaries it is an 
unfruitful task. The simplest way of being certain of the existence of a black hole with 
any u — > asymptotics is the following. 

First, let us recall that in w-coordinates the horizon is located at infinity. Therefore, 
we face the difficulty of estimating how representative some given numbers are of the 
actual solution near the horizon. Fortunately, assuming that the nature of the horizon 
is of the type of nonextremal D3-brane we can use formula ()2.32|) to estimate how close 
to the horizon we are. Numerically, we will consider a sequence for which p approaches 
the horizon as p = rjpo. Assuming that p$ = 2a we find that 

2u v a = -ln(l -rT 4 )- (3-8) 

This quantity characterizes how close we are to the horizon. Note that 

P^*= v -1. (3.9) 
Po 

Alternatively, we can look at this formula as giving a prescription that tells us for 
different values of a how far we have to go in u to reach the same "distance" from the 
horizon, where we refer to distance as a concept on p. Hence, we use expression (13. 8|) 
to estimate the value of u corresponding, for instance, to the 90% of the distance to 
the horizon. 

Next, we set the initial conditions for our variables at this value of W9o%, using 
expressions ()2.18|) for the solution of the non-extremal D3 case. We then integrate 
the system forward and identify numerically the presence of a horizon. To do that we 
developed the following set of criteria. 

Let us start by considering a massive geodesic in the general metric (j2.1j) . The 
effective Lagrangian for radial motion is given by 

C = -e 2z ~ 6x i 2 + e~ 2z+Wy u 2 . (3.10) 

Normalizing the effective Lagrangian to be £ = — 1 and using the fact that the La- 
grangian does not explicitly depends on time we find an equation for the radial coor- 
dinate: 

u 2 + e 2z - 10y = E 2 e 6x - 10y . (3.11) 

This equation does not readily allows an interpretation in terms of classical motion of 
a particle with energy E in a potential. For this aim, we find it useful to define a new 
variable v given by dv = e~ Sx+5y du. The geodesic equation then becomes 

v 2 + e 2z ~ 6x = E 2 } (3.12) 
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where the functions x,y and z are now viewed as functions of v. In this presentation 
we can consider the function e 2z ~ 6x as an effective potential describing the classical 
motion of a particle. Using this analogy and that dv/du > 0, we conclude that if the 
potential has a wall at some large value of u it indicates a singularity and potential 
repulsive behavior of the form characteristic to the negative mass Schwarzschild or 
KT. It thus works as a necessary condition for the existence of a black hole that the 
potential develops no wall for large values of u. In fact, the case of the regular non 
extremal D3 brane shows that the potential actually vanishes near the horizon. We will 
thus, use the vanishing of the effective potential as a first signal for a possible horizon. 

Vanishing of the effective potential is not sufficient to declare the existence of a 
horizon at that point. We have to take into account the fact that basically goo ~ 
e 2z-6x an( j therefore the potential can go to zero and the asymptotic time needs not 
diverge. The possibility is given by some integrable singularity in the expression for 
the asymptotic time. 

This way, to clarify the existence of a horizon we next consider the proper and 
asymptotic times. For a massive particle with energy E we obtain that the affine 
parameter is given by 

t= / =du, (3.13) 

J Ul VE 2 e^ - e 2z V ' 

On the other hand, the asymptotic time is 

t= / e~ 2z+6x du. (3.14) 

J Ul VE 2 e 6x - e 2z 

Since we are assuming an AdS configuration, in the above definitions there are some 
subtleties regarding the choice of the lower limit to ensure that the particle is in the 
energetically allowed region. However, this will amount to a numerical correction which 
will not modify the convergence of the above integrals when u — ► 00. 

To have the asymptotic infinity causally disconnected from the region behind the 
horizon we demand the affine time to be finite while the asymptotic time must be 
unbounded. For the numerical calculation of the above integrals we decided to imple- 
ment a version of the 3-5 points adaptive Simpson's rule intended for controlling the 
accuracy while making an optimal use of the output procedures calling for the solution 
of the differential system. 

Finally, we check for the the finiteness of the horizon area. Note at this point that, 
according to the analysis of subsection 12.2.21 we have the possibility of doing a cross- 
check of the accuracy of the numerical output as well as of the assumptions made for 
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that analysis. What we have to do is to observe the asymptotics of variables y and 
z and, if the horizon area and the corresponding temperature are finite, then y and 
z must exhibit a linear behavior when u — > oo, and, for a given a, the values of the 
corresponding slopes should obey condition (|2.4(J|) . 

Having found a regular black hole, we integrate backward toward u = 0. In order 
to have an idea of the analytical behavior hidden behind the numbers, we try several 
different fits by solving least-squares problems of the difference between objective func- 
tions and the numerical results. In the case of variables with regular behavior at the 
origin, we find the coefficients of a truncated Taylor expansion. In those cases with 
singular behavior, we use a set of singular real functions which includes the logarithm, 
the square root and rational functions. As a cross check, we also looked for the best 
fit of the corresponding derivatives. 

4 Numerical solutions with varying flux and nonextremal D3- 
like horizon 

Following the recipe in the previous section we were able to find several solutions with 
varying flux and nonextremal horizon of the form of D3 branes. 

4.1 Behavior near the horizon 

Recall that we are using the nonextremal D3 solution to fix the boundary conditions 
at ugo%. In other words, first we look for a certain point of the state-space (the 
space spanned by our ten variables) crossed by a trajectory when P = 0. Next, 
we assume that, after setting P / 0, close enough to that point there still remain 
trajectories leading to a black hole solution. We guess the position of a new point in 
one of these trajectories by looking to the corrections of the original point coordinates 
necessary for the Hamiltonian constrain (j2.10j) still being satisfied. It turned out that 
this methodology is very effective on the search for regular black holes. 

In Figs. El IU El an d El solid black curves represent the metrics functions obtained 
numerically for P = a = 1000 and Q = 1. We have chosen these values of the 
parameters to make sure the solution is quite different from the D3 or KT solutions. 
In these same figures, the corresponding nonextremal D3 solution, i.e., that with P = 0, 
have been drawn with red dashed curves. As it can be observed, the numerical solution 
differs from the D3 one. The difference is more pronounced for u — > 0, but in Figs. 
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Figure 3: Plot of the goo component of the metrics Ij'i.lj) . The solid black curve repre- 
sents the numerical solution with P = a = 1000 and Q — 1, and the dashed red curve 
the analytical solution with P = 0, a = 1000 and Q = 1. Here m 90 % ~ 0.0001. 
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Figure 4: Plot of the gu component of the metrics (|2.1|) . The solid black curve repre- 
sents the numerical solution with P = a = 1000 and Q = 1, and the dashed red curve 
the analytical solution with P = 0, a = 1000 and Q = 1. Here u 90 % ~ 0.0001. 
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Figure 5: Plot of the g uu component of the metrics 1)2.1)1 . The solid black curve 
represents the numerical solution with P = a = 1000 and Q — 1, and the dashed red 
curve the analytical solution with P = 0, a = 1000 and Q = 1. Here u$q% ~ 0.0001. 
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Figure 6: Plot of the #M5 component of the metrics (j2.1j) . The solid black curve 
represents the numerical solution with P = a = 1000 and Q = 1, and the dashed red 
curve the analytical solution with P = 0, a = 1000 and Q — 1. Here Uqq% ~ 0.0001 
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Figure 7: Plot of the logarithm of the effective potential. Here Ug % ~ 0.0001. 

E] and El we note that even close to the horizon there is some significant divergence 
between both solutions. Here we need to take into account that for these values of the 
parameters m 90 % ~ 0.0001. 

Nevertheless, the numerical results certainly indicate the existence of a black hole. 
To show that, in Figs. l7]l8|l9landlTUl we present, respectively, the natural logarithm of 
the effective potential Ve// = e 2z ~ 6x entering equation ()3.12|) . the proper and asymp- 
totic times as given by equations ()3.13|) and ()3.14|) . and quantity e - 2z + 3x + 5 v which, as 
it was mentioned in section I2.2.2[ is proportional to the black hole horizon area when 
u — > oo. 

Indeed, the following key properties are observed: absence of a potential barrier in 
the effective potential, the affine parameter converges to a finite value while the test 
particle approaches the horizon, simultaneously the asymptotic time diverges, and the 
corresponding horizon area is finite. Thus, we conclude that this solution represents a 
regular black hole for u — > oo. 

As discussed earlier in subsection 13.11 we can assess the quality of these results by 
analyzing the asymptotic behavior of fields y and z, which are plotted in figure ^2 
It is observed that both fields behave as linear functions of u. The linear behavior is 
required from the existence and regularity of the horizon, see discussion around (|2.39|) . 
As a matter of fact, a least-squares fit to a polynomial of this numerical output yields, 



2/oo = 1.90027 - 999.942947m + 0.0000002 u 2 + 0.00000002 u 3 + 0.000000002 u A , (4.1) 
Zoo = -0.000113 - 999.857345 m + 0.0000004 m 2 + 0.0000001m 3 + 0.00000001 w 4 . (4.2) 
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Figure 8: Plot of the affirie parameter. Here m 90 % ~ 0.0001. 
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Figure 10: Plot of exp (— 2z + 3a; + by). In the limit u — > oo this quantity is 
tional to the horizon area. Here u 90 % ~ 0.0001. 
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Figure 11: Plot of the fields y (black) and z (red). Here u 90 % ~ 0.0001. 
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We see that higher order coefficients are negligible. Using that here a = 1000, we find 
that a = -0.999857 and (3 = -0.999943, yielding -2a + 5(3 = -3.00000004. This is 
indeed, a strong result supporting the reliability of the numerical solution, as well as 
the corresponding asymptotical analysis presented in the last part of subsection 12.2.21 
According with expression (|2.44j) and the above results, the temperature of this 
black hole is equal to 0.035695 in the appropriate units. 



4.2 Asymptotical behavior at infinity, u = 

While studying the asymptotic behavior of the numerical solution at infinity, u — > 0, 
we have found that typically the z-field diverges at some u = u S i ng , while the remaining 
fields seem to be analytical at that point. Numerically, the value of u S i ng can be positive 
or negative. If u sing < 0, it means that the corresponding space-time is complete, and 
it is asymptotically flat. If u s i ng > 0, the corresponding space-time is singular in the 
sense that it collapses at a finite value of the radial coordinate. In principle, there 
exists also the possibility of u s i ng = 0. In this last case, the space-time is complete 
and it has a potential wall at infinity, like in AdS. Independently of the value of u S i ng , 
the qualitative behavior of the solutions is essentially the same. As we will see soon, 
obtaining one case or the other depends mainly on the value of a, as well as, on the 
boundary conditions. 

The interesting result here is that this singularity has a lot in common with the one 
in the KT solution. If we assume that, 



sinq 1 i 

— In 



Au In 



u 



Using 



4 

where A is a constant, then the corresponding derivative 

rl-y 1 l n 



(4.3) 



du 4 uln 



(4.4) 



depends only on u S i ng , which is given by the numerical method. In figure ^] the 
diamonds represent the numerical solution for dz/du with P — 1, Q — 1, a — 0.01. For 
this case u 90 % ~ 173 and we found u S i ng ~ 4.78571. As it can be noted, the behavior 
behind the numbers is very accurately reproduced by the solid curve corresponding to 
the plot of the function given by equation (|4.4jl . The corresponding behavior for z(u) is 
shown in figure El The main difference between these expressions and the KT solution 
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Figure 12: Understanding the singular behavior near u = u S i ng . The diamonds repre- 
sent the numerical output of dz/du, while the curve is the plot of the function given 
by equation ([4.4)1 . Here ugo% ~ 173 and u S i ng ~ 4.78571. 
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Figure 13: Understanding the singular behavior near u = u sing . The diamonds rep- 
resent the numerical output of z, while the curve is the plot of the function given by 
equation IjOl . Here u 90% w 173, u sing w 4.78571 and A w 1.6. 
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is that in the latter the range of u is < u < u S i ng , while here oo > u > u S i ng . This 
requires a change in the sign of A and, indeed, we have found that here A ps 1.6 > 0, 
while for KT, A = —P 2 /2. To confirm these findings, we fitted the solutions for y and 
/ (and their derivatives), using respectively the following expressions, 

1 2 

y ps -7 In \Au\ + y.O + y_i (u - u sing ) + y_2 (it - u sing ) , (4.5) 

+ f_0 + f-1 (u - u sing ) + fJ2 (u - u smg ) 2 , (4.6) 

where e is the Euler number. Solving a least-squares problem for the coefficients, we ob- 
tained y_0 = -0.041239, yA = 0.07182457, y_2 = -0.128777, B = 0.905268496, f_0 = 
1.274356, fA = 0.285765, f_2 = -0.05121489. The fits are amazingly good, with 
X 2 < 10 -10 for both cases. Since, these variables are analytical at u S i ng , we tried a 
direct fit with Taylor expansions for them, as well as for the remaining fields. The 
solution of the corresponding least-squares problems yielded, 

y(u) ~ -0.783602 - 0.042953(m - u smg ) + 0.003493(m - u sing ) 2 

+ 0.001668(m - u sing f - 0.115651(n - u smg ) 4 + 3.027465(u - u sing f ',(4.7) 

f(u) ~ 0.34785 - 0.061501(u - u sing ) + 0.008764(u - u sing f 

- 0.002828(w - u sin gf + 0.099713(m - u smg f - 2.972385(u - u sing f ,(4.8) 
$(w) ps 0.231423 - 0.027829(w - u sing ) + 0.002474(w - u sing f 

- 0.002109(m - u S i ng f + 0.120818(m - u S in 9 f - 3.620755(m - M sin5 ) 5 ,(4.9) 
w(u) ps 0.028572 - 0.004287(m - w sm9 ) + 0.000513(m - u sing f 

0.000065( ? i - u smg f + 0.000009( M - u sina f - 10~ 6 (m - u sing f. (4.10) 

We recall that these results were checked by also fitting the corresponding derivatives. 
Expanding the non-polynomial part of the expressions (|4.5J1 and (J4.6)) it is not difficult 
to test the direct correspondence between both approximations of the behavior of y 
and / when u u S i ng - 

Substituting in system ()2.9|) the ansatz ()4.3|) for the z-field and Taylor-like parametriza- 
tions for the remaining fields near u S i ng , it can be proved that they provide an asymp- 
totic solution if the following condition is satisfied, 

f — -J_p e ($0+4yo+4wo) 

We can check if this analytic result is verified numerically. Indeed, using the numbers 
in the Taylor expansion of the fields given above from the numerical analysis, and 
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P = 1, we can verify that 



0.061501 « /i 



e 0.231423+4(-0.783602)+4(0.028572) 




Thus the analytic condition is indeed fulfilled by the numerical analysis. As a matter 
of fact, we have found that this condition is satisfied by all the solutions we analyzed. 
For instance, the solution with P = Q = 1 and a = 0.1 is regular at the origin, because 
Using — —0.056996 which corresponds to an asymptotically flat solution . For this case 
we also verify that, 



This is a very strong validation of our numerical results as well as of the analytical 
solution for the u — > u sing asymptotics given by expressions (jOJ), f)4.5|) and ([4.6)1 . 

From the point of view of field theory the more interesting solution is the one with 
u S ing = 0. We have found that the position of the singularity near the origin depends 
mainly on the values of a and of the boundary conditions. The dependence on a seems 
to be described by a convex function u S i ng (a) having a minimum at some a = a c . 
This minimum will be positive or negative depending on the boundary conditions. 
There is a unique set of boundary conditions that lead to u S i ng (a c ) = 0. The main 
difficulty to find such a solution is that, because of the finitness of the integration step, 
falling into a point with u S i n g < 10~ 10 is very difficult when the integration starts far 
from that value. We should also recall that ours is a nonlinear system, therefore it is 
sensitive to changes of the parameters and strongly sensitive to changes of the boundary 
conditions. There is also a restriction on varying the boundary conditions coming from 
the necessity of preserving the black hole solution at u — > oo. Nevertheless, we were 
able to obtain solutions near enough to the origin (for instance, with \u S i ng \ = O(10 -6 )) 
and independently of the sign we have always found that the whole analysis in this 
section applies, so there is no particular reason to expect the situation to be different 

for Using = 0. 

To finish this section we would like to remark that we observed that after reaching 
the minimum, u S i ng gets larger while still increasing a. On the other hand, u 90 % always 
decreases while increasing a. So, there exists a not necessarily high value of a such 
that Using = ugo%. As a matter of fact, that is true for any distance until the horizon as 
estimated by expresion (|3.8|) . It implies that there seems to always exist a maximum 
value of a such that there is no solution outside the black hole horizon. 



-37.008 « / : 



e 0.36724+4(0.693744)+4(0.117277) 



37.015. 
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4.3 Area and temperature dependence on the parameters 

In this paper we will not attempt a full analysis of the dependence of the solution on the 
three parameters a, P and Q, as well as on the boundary conditions. We concentrate 
on explicitly constructing numerical solutions. Nevertheless, we would like to present 
some observations that show the power of the numerical method to understand the 
physics behind the more general solutions. The naive intuition is that the properties 
of the solution are controlled by parameters like 

P 2 a a , 1 \ 

~Q' Q- W (4U) 

This was however, not supported by the numerical analysis. 

As we already mentioned fixing P = always yields the nonextremal D3 with 
regular horizon. Increasing P from zero the solutions differs very slowly from the D3 
solution. For values lower than P = 10~ 6 , the existing black hole solutions cannot be 
practically distinguished from the D3 one. Nevertheless, it must be noted that some 
intervals of P arise where black hole solutions cease to exist. For instance, we failed 
to find a black hole for Q = and a = O.l and P = 10~ 8 . Above some value 

of P, the horizon area does not seem to converge to a finite value. And this is not 
necessarily a very high value of P. For instance, for Q = 10~ 4 and a = O.l, as above, 
and P > 5 x 10~ 2 , we found nothing but degenerate black hole solutions. So, we are left 
here with the range P G [10~ 6 ,5 x 10~ 2 ), where the dependence on this parameter is 
very weak. In principle, as one increases the value of P the temperature also increases. 
The horizon area has the same behavior. On the other asymptotic regime of u, for 
small u, the larger P, the lower the value of u S i ng seems to be. However, this behavior 
is altered near those islands where solutions with regular horizons do not exist. 

The dependence on Q appears to be the opposite of what we observed for P. For 
instance, with P = 10~ 7 and a = O.l, for Q < 10~ 10 the horizon area does not converge 
to a finite value. Increasing Q up to 4.4 x 10~ 4 some intervals arise where black hole 
solutions exist and others where they do not. If a black hole solution exists, then for 
all value of Q it has a finite horizon area. For 4.4 x 10~ 4 < Q < 10 20 all the solutions 
we found were regular black holes. As in the D3 solution, Q does not seem to affect 
the asymptotic behavior of y when u goes to infinity, but affects the asymptotics of 
z. As it is shown in figure El Q is inversely proportional to the temperature. The 
horizon area has the inverse behavior. This kind of dependence are proper also of the 
D3 solution. The value of u sing does not seem to depend on Q. 
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Figure 14: The dependence of the black hole temperature on Q. The shaded regions 
mark some intervals of Q where no black hole solution seems to exist. 

The first thing we note while varying a is that, given a combination of P and Q, 
no value of a exist (except a = 0) that could change one of the following results: 
no black hole, singular black hole, regular black hole. The larger a, the higher the 
temperature, which is also the kind of dependence observed in the D3 solution. In 
figure ITB1 we present the result for P = 10~ 7 and Q — 1. The horizon area presents a 
similar behavior. Figure suggests that the horizon temperature and area blow up 
at some finite value of a. That confirms the findings described at the end of section 
14. 2\ where it was discussed the dependence on a of the critical value u S i ng . It means 
that there seems to exists a maximal value of a above which no black hole with finite 
temperature or area can exist. 

5 Conclusions 

In this paper we have solved numerically the IIB supergravity equations of motion in 
the presence of varying flux. This class of solutions is controlled by three parameters 
loosely identified as Q with the five-form flux, P with the three-form flux and a with 
nonextremality. 

We established that choosing the right boundary conditions, and for a wide range 
of values of the parameters P, Q and a, we can obtain solutions representing black 
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Figure 15: The dependence of the black hole temperature on a. 



holes with finite horizon area. The criteria to assert the existence of the black hole are: 
no barrier of the effective gravitational potential, convergence of the affine parameter, 
divergence of the asymptotic time and finiteness of the horizon area. In this regime 
the obtained solutions resemble but are not identical to the standard non-extremal D3 
solution. 

We also detected some curious behavior which deserves a more systematic study. 
Let us summarize how the horizon temperature and its area depend on the parameters 
of the solution. In the parameter space (P, Q) there are islands of values which do 
not lead to black hole solution. This particular result weakly depends on the values 
of a. When a black hole exist, the temperature and area of the horizon seem to be 
proportional to the value of P. Above some value of P, the horizon area ceases to 
converge to a finite value. The dependence of the horizon temperature and the area on 
Q is very much like that in the D3 solution, i.e, the higher Q, the higher the horizon area 
but the lower its temperature. Curiously, below some threshold value of Q the horizon 
area is no longer finite. Finally, the dependence of the thermodynamical quantities 
on a is quite interesting. It resembles the situation in the D3 solution, the area and 
the temperature are both proportional to a. However, here our results indicate the 
possibility of the existence of a maximal value of a above which finite-temperature 
black hole can not exist. 

In the region asymptotically away from the horizon we were able to obtain a very 
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accurate description based on a combination of analytic and numerical methods. We 
found that the solution has qualitatively a KT-like asymptotic behavior but with mod- 
ification allowing the solutions to exist for u > u sing . Depending on whether u sing is 
negative or positive we have an asymptotically flat or a collapsing space-time. There 
seems to exist a solution with u S i ng = which is complete and behaves asymptotically 
as the AdS space-time. 

Let us end our remarks about the solution by presenting the dependence of u sing 
on the parameters of the solution. While increasing a, u S i ng initially decreases, reaches 
a minimal value and then starts to increase. This implies the existence of a maximal 
value of a above which no black hole with finite temperature or area can exist. The 
dependence of u S i ng on the remaining parameters is more complex and deserves further 
study. Though weakly, the value of u sing seems to be inversely proportional to P and 
almost independent of Q. 

In this paper we have focused on constructing the solution. There are many ideas 
that would be interesting to explore on this supergravity background. An important 
venue that we plan to return to [21] is the possibility of a Hawking-Page phase transition 
between the solution constructed here and the Klebanov-Strassler background. Some 
of the various questions we plan to address in future publications have been addressed 
recently in simpler supergravity backgrounds |2"2"1 1231 
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A Geodesic analysis of some typical solutions 

In this appendix we explicitly discuss the geodesic structure of some known solutions. 
The idea is to develop the necessary intuition to interpret the numerical results. We 
start by considering geodesic motion in the general metric (j2.1j) . The effective potential 
for radial motion is given by 

C = - e 2z - 6 H 2 + e- 2z+10y u 2 . (A.l) 
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Normalizing the effective Lagrangian to be C = — 1 and using the fact that the La- 
grangian does not explicitly depends on time we find an equation for the radial coor- 
dinate: 

ii 2 + e 2z - l0y = E 2 e 6x - 1Q y. (A.2) 

This equation does not readily allows an interpretation in terms of classical motion of 
a particle with energy E in a potential. For this aim, we find it useful to define a new 
variable v given by dv = e~ 3x+5y du. The geodesic equation then becomes 

v 2 + e 2z - 6x = E 2 , (A.3) 

where the functions x, y and z are now viewed as functions of v. In this presentation 
we can consider the function e 2z ~ 6x as an effective gravitational potential describing 
the classical motion of a particle. 



A.l AdS space 

For AdS space in Poincare coordinates which are the natural coordinate for the D3 
brane we find 

The equation for radial motion as a function of the proper time r is: 

r 2 + ^ = E 2 . (A.5) 

The corresponding potential is V(r) = r 2 /L 2 . This equation can be easily recognized 
as the harmonic oscillator whose solution is 

r(r) = EL sin(y r + 0), (A.6) 

where is a phase that determines the initial condition for the motion of the particle. 

Our main message is that "physical" experiments in AdS are quite different from 
Schwarzschild. As we can see from the potential, AdS functions acts as a box and the 
experiments we can perform involve sending a particle from the interior of AdS towards 
the boundary r = oo. Such particle explores the boundary and returns back having 
reached its maximum radius r max = EL. 
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A. 2 Nonextremal D3 branes 

For the nonextremal D3 brane in the standard parametrization the equation for massive 
geodesic is 



2 4 

f'+U-it^V- (A.7) 
Note that taking r = yields an equation which we recognize as the harmonic oscillator 
describing motion of massive particles in AdS. More generally, we can view this problem 
as a classical mechanical problem in a potential given by 

2 4 

For large r we have the typical harmonic potential of AdS but for r — > r the potential 
goes to zero. 

The "physical" experiment here is slightly different. We have motion between r = tq 
and some r max given by: 

,4 



(EL) 2 = r 2 max + (A.9) 



„2 
' max 

Note that r max in this case is larger than EL which corresponds to AdS with no black 
hole. Namely, 

The massive geodesic in the p coordinates has 

dp 



T, 



massive 



dp EL 2 



L 2 



A. 3 The structure of the KT solution 



£ 2 _|i !_4 H 1 



.*\ *.2 M _ £o 



4\ ' 



(All) 



Let us verify that the KT solution has a behavior similar to the negative mass Schwarzschild 
solution and that it can be detected by the effective potential (jA.3|) . Consider a solution 
of KT type with metric given by 

ds 2 = h- l l 2 dt 2 + h l l 2 dr 2 + ... (A.12) 
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where the warp factor is of the form 



h=*{l + P\n-), (A.13) 

r 4 r 



A massive geodesic follows an equation of the form: 



2 

r 2 + . r = E\ (A.14) 

L\ 1 + Pln^ 



ro 

Let us, once again, consider the natural limits. For large values of r the potential is 
mainly harmonic V(r) ~ r 2 /li)V 2 (r). This behavior is typical of asymptotically AdS 
spaces as we saw at the beginning of the section. For the purpose of a horizon we are 
interested in whether or not the potential develops a wall for small values of r. Indeed, 
it is clear that the denominator of the potential is zero, that is, near 

r S ing = r exp(-l/P). (A. 15) 

Around this point the potential goes to infinity and the spacetime behaves like in the 
neighborhood of a negative-mass Schwarzschild black hole. Before this point there is 
a minimum of the potential located at 

P — 4 

r m in = r exp( — ) = r sing e 1/4 . (A. 16) 
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